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It has been proposed that multifragmentation can be related to the 
hqiiid-gas phase transition of nuclear matter. We study the statistical 
properties of finite nuclear matter near the phase transition with the 
help of a Lattice Gas Model (LGM) . The original version of LGM with 
only one type of charge-neutral particles is well known to feature the 
properties of the liquid-gas phase transition. In this contribution, we 
address the effect of Coulomb and isospin dependence interaction for 
the finite nuclei transition, and study the symmetry energy properties 
of finite temperature systems. 

In the last decade, there has been a growing interest for the measurement 
of the symmetry coefficient Csym of the nuclear equation of state, through 
experiments with heavy ion coUisions. 

At zero temperature, most mean-field-based nuclear equations of state 
can be approximately reproduced through the density (p) functional form 



where pQ is the saturation density, and 7 is a coefficient which deter- 
mines the stiffness of the symmetry coefficient. The values of 7 can widely 
differ among different EOS parametrisations (0.5 ^ 2) and are not presently 
constraint by experimental measurements. 
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One proposed technique is the so-called isoscaling analysis. It is well 
known [2] that the ratio R2i{N,Z) of isotope yields Yi{N,Z) measured in 
two reactions similar in temperature but different in isospin labelled (1,2) 
has an exponential behavior according to : 



where a and /? are the isoscaling parameters. 

Different data have been published showing that a decreases with in- 
creasing the violence of the collision [21 |3l |3] . Since an increasing collision 
violence is associated to increasing temperature and/or decreasing density 
at the time of fragment formation, this suggests that we may access the den- 
sity and/or temperature evolution of the symmetry term. Indeed at thermal 
equilibration and in the framework of the grand canonical ensemble, a and 
Csym are approximately related by [5] : 



where T is the common temperature of the two systems, and Zi,Ai are 
the corresponding charges and masses. 

Heavy ion collisions are probably the best probe of nuclear properties in 
excitation energy and density conditions far from the ground state. How- 
ever, the connexion between the symmetry energy coefficient appearing in 
eqsll] and [3] is far from being clear. Indeed the fragmented configurations 
accessed by isoscaling analysis are dishomogeneous and present high order 
correlations; it is therefore not clear whether, even at thermal equilibrium, 
the associated energy functional only depends on the global density and tem- 
perature as obtained in a mean-field based picture. Moreover eql3]has been 
derived in the framework of macroscopic statistical models l5j where many 
body correlations are supposed to be entirely exhausted by clusterisation 
and clusters are described as independent degrees of freedom. If the validity 
of isoscaling eql2]in microscopic theories is well settled [6], the connexion of 
a and Cgym eqlU in such models has never been proved. 

To progress on these issues, it is interesting to consider a microscopic 
model simple enough to be exactly solvable through Monte-Carlo simulations 
without any mean-field or independent cluster approximation. 

In this contribution, we propose to study the temperature and / or density 
dependence of the symmetry coefficient Csym in a lattice gas model. 
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1 Description of the model 

We use a cubic lattice of linear site L = 20, where each cell i is characterised 
by four degrees of freedom : one discrete variable ai for isospin (fip = 1 , 0"^ = 
— 1, ctq = 0) and three continuous variables pi for the momentum. 
The hamiltonian of the system follows: 

(ij> a,=aj=l,{i^j) i=l 

where < i,j > are nearest neighbour cells, ecr-o-j is the coupling be- 
tween nearest neighbour (ei = 0,e_i = h.bMeV), Ic{= lAAMeV / fm) is 
the Coulomb interaction between all protons in the lattice, and rij is the 
distance between sites i and j. The last part of the interaction is the kinetic 
term where m = 939Mey is the nucleon mass. 

Calculations are made in the isobar canonical ensemble, which has been 
shown to be the correct canonical ensemble to describe unbound systems in 
the vacuum [7j. 

The partition fonction is then: 

(n) 

where the sum runs over all the possible realisations of the system, and 
^3(n) jg ^]^g global extension of the system for each partition (n) define as: 

The statistical average at a given value of T and P ^ill be 

noted in the following. 



2 Phase transition 



In its original version with only one type of particles and closest neighbours 
interaction, the Lattice Gas Model (LGM) is well known to exhibit a first 
order transition line at low temperature ending to a second order critical 
point [H]. In this case, with a short range isovector coupling and a long 
range repulsive interaction, the phase diagram is not qualitatively modi- 
fied [To]. Indeed we find a phase coexistence line ended by the critical point, 



characterized by its temperature and pressure Pc- In the following, all 
calculations are made at the subcritical pressure P = 5. x 10~^MeV.fm~^. 

In table [H we compile the transition temperature Tt for different systems 
at the chosen pressure. The inclusion of Coulomb and isospin dependence 
results in a lower transition temperature than for a scalar interaction without 
coulomb. With the realistic hamiltonian, in all cases the isospin dependence 
is weak on the transition temperature, less than few percents. 



System 


(75, 75)sc,alar 


(75,75) 


(85,65) 


(88,62) 


(91,59) 


TtiMeV) 


3.17 


2.27 


2.285 


2.290 


2.295 



Table 1: Liquid-Gas transition temperature of different systems oi A — 150 nuclcons 
at subcritical pressure P — 5.10~^ MeV.fm~^. 

To explore the finite temperature symmetry energy of the model, we will 
now use a macroscopic parametrisation in order to connect the microscopic 
and macroscopic properties of these systems. 

3 From microscopic to macroscopic 

As shown in the previous section, finite temperature Lattice Gas systems 
are strongly dishomogeneous and clusterized because of the presence of the 
phase transition. In this situation, it is not clear whether the global ener- 
getics of the system, including its symmetry properties, can be described by 
a macroscopic parametrisation depending on the average density as in the 
nuclear mean- field. 

To explore this issue and in order to extract the symmetry energy of 
our systems, we use a liquid-drop (macroscopic) parametrisation for the 
interaction energy of the system, which reads: 

El;^'^{6,p,T) = a,{p,T)A + as{p,T)A^/^ + Csym{p,T)A6^ + ac{p,T)^ (7) 

Here 6 is the isospin asymmetry (A^ — Z)/{N + Z), T is the temperature 
and p = A/ {A/ 2,71 R^) is an estimation of the average density of the system, 
where the mean cubic radius from eqi6]is calculated, excluding the monomers 
(^ = 1). 

This parametrisation uses four macroscopic parameters: is associated 
to the volume energy, corresponds to the surface energy, Cgym is the 
coefficient related to the symmetry energy, Uc corresponds to the Coulomb 
contribution of this interaction. 



3.1 Temperature dependence 



In order to extract all these parameters, we fit our systems at different 
temperatures. The results obtained from the parametrisation are plotted 
versus the true values coming out from the simulation on the figure [T]-/e/t 
We observe a good agreement between the two results (close to the line 
= Eint)- The values of the macroscopic coefficients are plotted as 
a function of the temperature on the figure fright. Oc is constant with 
temperature, while and Ug decrease with increasing temperature, and c^yrn 
decreases with increasing temperature with a bump around the transition 
temperature for the considered pressure (Tt ~ 2.25MeV). 




Figure 1: Left: Correlation between the liquid drop parametrisation (eq|7]) and the 
average LGM energy for eight different systems at different temperatures. Right: 
Evolution of the macroscopic coefficients with the temperature. 

A very strong temperature dependence is observed, at variance with 
nuclear mean-field calculations [1]. However, up to now, we have extracted 
the temperature dependence of the macroscopic coefficient without looking 
at the density dependence, while a temperature variation induces a change 
in the average system size for these constant pressure calculations [llj . 

3.2 Density dependence 

On figure ^left, we plot the evolution of the density with the temperature. 
We observe a decrease of the density with increasing temperature similar to 
the behavior of and a^. 

The temperature variation of and can be easily understood as a 
simple effect of this density change, as we now explain. 

To avoid the interference of the coulomb and symmetry effects, we con- 
sider an isospin symmetric system (N = 75, Z = 75) and switch off the 



coulomb interaction. 

To distinguish the respective role of p and T, we sort out the temperature 
dependence induced by the density variation p{T) as: 

avAP^ T) oc /^,s(p(T))a; ,(/9, T) (8) 

If the temperature dependence shown in figure[T]is only due to the change 
with temperature of the average size as the mean-field approximation, eqlH] 
should be fulfilled with a'^ and a'g as constants. This is indeed what is 
observed once the total bulk energy is plotted as a function of the mean 
radius, as shown in figure fright. The observed functional dependence of 
f{p) can be easily understood from simple geometrical considerations : 



/.(p(T)) a oc /^(p(T))cx|g) oci2-i(T) (9) 

where S and V represent the average system surface and volume. 
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Figure 2: Left: density evolution in fonction of the temperature. Right: Total bulk 
energy in fonction of the radius of the system. 

A similar analysis for the symmetry coefficient with the realistic hamil- 
tonian including Coulomb is currently under progress. 

Conclusions 

In order to connect the microscopic properties of the nuclear equation of 
state (symmetry energy and its dependence on temperature and/or density) 
and experimental fragment observables like isoscaling, we have used a lattice 
gas model with an isospin dependent short range and a long range coulomb 
interaction. In this case, we observe a liquid-gas transition at lower tempera- 
ture than the liquid gas transition with a purely nuclear isospin independent 
interaction. The coupling of an isovector and coulomb interaction induces 



another first order transition, a kind of fission-fusion transition at lower tem- 
perature than the Hquid-gas one. We have also established a macroscopic 
parametrisation of the model, with four parameters which are function of 
temperature and density (a 

V 1 (^s iCsymi Q^c )• For isospin symmetric systems, 
and as long as the coulomb interaction is ignored, the whole temperature de- 
pendence appears to be exhausted by the total average density dependence, 
even in the phase transition region where the system is highly dishomo- 
geneous and clusterised. This finding suggests that it may be possible to 
access the density dependence of Cgym in multifragmentation experiments, 
as it is actively searched for in several recent experimental analysis O HJ [5] . 
The connection between isoscaling observables and the symmetry energy for 
finite temperature within this model is currently in progress [iQi . 
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